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ABSTRACT 

We fir st show that the new WMAP 3 year data confirm the detection bv lMvers et alj 
(2004) of an extended SZ signal centred on 606 Abell (ACO) clusters with richness 
class, R ^ 2. Our results also show SZ decrements around APM and 2MASS groups 
at increased sign ificance than previously detected. We then follow the approach of 
iLieu et al . (2006) and compare the stacked WMAP results for the decrement in 31 
clusters with ROSAT X-ray profiles where Lieu et al found on average less SZ decre- 
ment in the WMAP 1 year data than predicted. We confirm that in the 3 year data 
these same clusters again show less SZ decrement than the X-ray data predicts. We 
then analysed the WMAP results fo r the 38 X-ray clu s ters w ith OVRO/BIMA mea- 
sured SZ decrements as presented by iBonamente et ail (|2006l ). We again find that the 
average decrement is measured to be significantly less (5.5a) than predicted by the 
Chandra X- ray data. Thus wh ile we confirm the original detection of an extended 
SZ effect by iMvers et al.l ((2004), these X-ray comparisons may now suggest that the 
central SZ amplitudes detected by WMAP may actually be lower than expected. One 
possible explanation is that there is contamination of the WMAP SZ signal by radio 
sources in the clusters but we argue that this appears implausible. We then consider the 
possibility that the SZ decrement has been lensed away by foreground galaxy groups. 
Such a model predicts that the SZ decrement should depend on cluster redshift. A 
reduction in the SZ decrement with redshift is suggested from the ACO cluster sam- 
ple an d also from comparing the samples of ILieu et ail (|2006h and IBonamente et all 
(20061). However, the mass power-spectrum would require a far higher amplitude than 
currently expected if lensing was to explain the SZ deficit in high redshift clusters. 
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1 INTRODUCTION 

I Myers et alj (|2004T ) made a cross-correlation analysis be- 
tween galaxy cluster catalo gues and the WMAP first year 
data ('Hi nshaw et al.l [20031 ) . They saw a statistical decre- 
ment near groups and clusters as detected by APM and 
also in more nearby groups and clusters as detected by 
2MASS but the strongest signal was seen in the ACO rich 
cluster catalogue. There the decrement was approximately 
what was expected from predictions based on X-ray obser- 
vations of the Coma cluster which is itself a richness class 2 
cluster. However, the profile appeared to be more extended 
than expected from simple fits to these typical cluster X-ray 
dat a. The exten t of th e SZ effect, possibly to 8 « 1 degree, 
led IMvers et ail (|2004l ) to speculate whether the SZ effect 
could contaminate the measurement of the acoustic peaks, 
although the difference between the SZ and primordial CMB 
spectral indices may constrain such a possibility at least for 
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the first peak (|Huffenberger et alJliooi ). We now return to 
this topic with the first aim to see if the extended SZ effect 
reproduces in the 3-year WMAP data. 

Meanwhile. ILieu et al.l (|2006l ) analysed the WMAP first 
year data now focusing only on 31 clusters with ROSAT X- 
ray data. They made basic predictions for the SZ decre- 
ment in each cluster and found that they over-predicted 
the SZ decrement. One possibility was that discrete ra- 
dio sources in the clusters w ere diluting the decrements 
but this was argued ag ainst bv lLieu et alj ([2006). However, 
ILieu fc Quenbvl (|2006h suggested an alternative mechanism 
based on synchrotron radiation from cosmic ray electrons 
moving in the cluster magnetic field forming a diffuse clus- 
ter radio source which again may dilute the SZ effect. This 
model was also aimed at explaining the soft X-ray excesses 
detected in some clusters via inverse Compton scattering of 
the CMB by the same cosm ic ray electrons in the cluster 
(e.g. iNevalainen et alj|2003l and references therein). 

Here we shall check the result of ILieu et al.l l|2006l ) us- 
ing our cross-correlation methodology and the full WMAP 
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Table 1. Properties of the WMAP frequency bands. 



Band 


Frequency 


FWHM 


W 


94GHz 


12/6 


V 


61GHz 


19/8 


Q 


41GHz 


29/4 


Ka 


33GHz 


37/2 


K 


23GHz 


49/2 



3-year da ta. In the f irst in stance, we shall take the X-ray 
models of lLieu et al.l l|2006h which follow the simple (3 model 
prescription described in Section 3 below. We shall also look 
at a new sample of clusters with excellent Chandra X-ray 
data (jBonamente et all 12006^ . Again we shall simply take 
their models convolved for the WMAP PSF in the appro- 
priate band and compare to the averaged SZ decrement seen 
in the WMAP3 data. 



The sky density of groups and clusters is 3.5 deg 2 at m ^ 7 
and 0.35 deg -2 at m ^ 15. We assume an average redshift 
of z = 0.1 for both APM samples. 

2.2.3 2MASS 

The third cluster catalogue is derived from the final data 
release of the 2M A SS Extended Source Catalogue (XSC) 
jjarrett et all l200d ) to a limit of K s < 13.7. ^-selected 
galaxy samples are dominated by early-type galaxies which 
are the most common galaxy-type found in rich galaxy 
clusters. Therefore the 2MASS survey provides an excel- 
lent tracer of the high density parts of the Universe out 
to z < 0.15 and so provides a further test for the existence 
of the SZ effect. Using the above 2-D friends-of-friends algo- 
rithm, I Myers et al.l l|2004 ) detected 500 groups and clusters 
with m ^ 35 members at the density contrast a = 8 in the 
|6| ^ 10 deg area. The 2MASS groups have average redshift, 
z « 0.06. 



2 DATA 

2.1 WMAP Third Year Data 

In this paper we use the raw CMB tempera ture maps pro- 
vided in the WMAP 3 year data release l|Hinshaw et al.l 
2006). These consist of temperature data from the five fre- 
quency bands and the internal linear combination (ILC) map 
(Table [TJ. In order to remove contamination from our own 
galaxy, we make use of the KpO foreground mask made avail- 
able with the other WMAP data products and have applied 
this to all our maps prior to cross-correlation. Finally, the 
data is used here in the HEALPix format of equal area data 
elements, characterised by N si[ j e =512, which gives an ele- 
ment width of « 7'. 

2.2 Cluster Data 

2.2.1 AGO 

The ACO catalogue (jAbell. Corwin fc Olowml 1 19891 ) lists 
clusters with 30 or more members, given the requirements 
that all members are within 2 magnitudes of the third 
brightest cluster member, whilst also lying within a 1.5 h _1 
Mpc radius. A richness class, R, is applied to the individual 
clusters based on a scale of ^ R ^ 5. The catalogue covers 
both hemispheres and here we trim these samples such that 
we take clusters of only R 2 and galactic latitudes of |b| 
^ 40°. 

2.2.2 APM 

We shall also use galaxy group and cluster catalogues de - 
rived from the APM Galaxy Survey of lMaddox et all (|l990h 
which covers the whole area with S < —2.5 deg and b < 
—40 deg. These were i dentified using the same 'friends- 
of-friends' algorithm as I Myers et al.l (|2003l ) and references 
therein. Circles around each APM galaxy with B < 20.5 are 
'grown' until the over-density, a, falls to a = 8 and those 
galaxies whose circles overlap are called groups. The APM 
galaxy surface density is N « 750 deg -2 at B < 20.5. Min- 
imum memberships, m, of m ^ 7 and m ^ 15 were used. 



2.2. 4 ROSAT X-ray cluster sample 

The 31 clusters published bv lBonamente et al.l |2002) were 
originally selected as a sample of X-ray bright clusters suit- 
able for observing X-ray surface brightness profiles. These 
profiles were obtained with the ROSAT PSPC instrument 
and estimates of the gas temperature, density and distri- 
bution were made by fitting a (3 profile model to the data 
(see section 3 below). T he X-ray data for these 31 clusters 
were previously used bv lLieu et al.l (|2006l ) to construct pre- 
dictive models of the SZ profile of each cluster. Redshifts for 
the clusters range from z~0.02 (Coma) up to z~0.3 (Abell 
2744), whilst the sample lies in the galactic latitude range 
of [ b |> 25°. 

2.2.5 Chandra X-ray cluster sample 

We further analyse the 38 clusters discussed by 
iBonamente et ail (|2006l ). These clusters h ave been observed 
at 30 GHz by OVRO and BIMA (see IBonamente et all 
(2006) and references therein) to detect the SZ decrements 
and have also been observed by Chandra to provide the 
X-ray data needed to estimate the value of Ho. The interfer- 
ometric radio observations have a resolution of ~ 1' and the 
X-ray observations from the Chandra ACIS-I camera have 
a resolution of ~ 1". Redshifts for these clusters are in the 
range 0. 18 < z < 0.8, a h ig her ra nge than for the ROSAT 
sample. IBonamente et ail l|2006h fitted both hydrostatic 
equilibrium and isothermal j3 models to the X-ray data and 
made predictions for the SZ decrements. 



3 SZ X-RAY MODELLING 

The SZ effect is generally modelled using X-ray gas profiles, 
densities and temperatures. The X-ray data is most simply 
modelled by fitting a /3 model to the X-ray intensity profile: 

/ a2\ (l-6/3)/2 

Sx = S X0 + ¥ J (1) 

where Sxo is the central X-ray surface brightness and 8 C is 
the angular core radius. On the isothermal assumption, the 
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Figure 1. Cross-correlation results between the WMAP 3-year W-band temperature data and the four cluster datasets: (a) 2MASS, 
(b) ACO, (c) APM m5;15, and (d) APM m^7. The dashed and dot-dashed lines in (b) show SZ models with AT = 0.083K" and 
AT = 0A9K respectively, both with 8 C = 1.'5 and j3 = 0.75 and convolved with the WMAP beam-width. The latter model is intended 
to be representative of the Coma cluster, scaled to rcdshift z = 0.15. The former is the ACO model fitted bv lMvers et al. 1 1120041) in their 
analysis of the WMAP 1st year results. 
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Figure 2. The cross-correlation of the ACO catalogue is shown 
after increments in galactic longitude of 20deg in the Abell cluster 
positions. The mean AT is shown for WMAP pixels within 6/3, 
60' and 500' of cluster centres, where the significance at each 
angular limit is 3.2(7, 2. Oct and 1.2ct respectively. 



temperature decrement, ATsz, as a function of the angular 
distance from the cluster-centre, 9, is then given by: 



AT SZ (6) = AT sz (0) 



1+ - 



.M. 



(2) 



Then the magnitude of the central temperature shift, 
ATsz(O), is given by: 



Ar sz (o) 

TcMB 



kT e 



dl 



x(e x + 1) 
e x - 1 



(3) 



where x = hvjkT^ , otk is the Thomson cross-section and 
n e , T e are the gas density and temperature derived from the 
X-r ay data. 



Lieu et all [2006) use the cluster sample of 



iBonamente et aT (2002) and fit ROSAT PSPC cluster X- 
ray profiles. Th ey assume isothermal ga s distributions with 
T e ta ken from IBonamente et all l|2002l ). IBonamente et all 
(2006) use both a hydrostatic equilibrium model, allowing 
a double power-law /3-model to allow for variations in the 
number density with radius, and an isothermal /3-model. 
With the hydrostatic model, they allow the gas temperature 
to vary with radius and a CDM component as well as gas 
to contribute to the cluster potential. We shall simply 



assume the isoth e rmal models of iLieu et ail (|2006l ) and 
IBonamente et all ([2006) and convolve the predicted SZ 
profile with the appropriate WMAP beam profile, modelled 
as a Gaussian with the FWHM beam-widths shown in 
Table □ 



4 CROSS-CORRELATION ANALYSIS 

We focus our analysis on the 94GHz W band from WMAP, 
looking for correlations characteristic of the SZ effect in this, 
the highest resolution b and. We perform a cross-correlation 
analysis as described in lMvers et al.l (|2004l ). calculating the 
mean temperature decrement/increment as a function of an- 
gular separation from galaxy clusters in the above datasets. 
Our cross-correlation takes the form: 



AT C (9) = J2 



ATi(fl) - AT 

Tli(0) 



(4) 



Where ATi(9) is the WMAP temperature in an element i 
at an angular separation 9 from a cluster centre and is 
the number of elements at that separation. AT is the mean 
WMAP temperature decrement across the entire region used 
in the analysis. For the 3-year W-band data, AT ~ 10~ 3 mK. 

Errors on our results are estimated using repeated 
Monte Carlo realisations of the cluster data. As each cat- 
alogue (i.e. 2MASS, ACO, APM) will each be highly clus- 
tered, it is important to incorporate this clustering into the 
realisations. Thus for each realisation we begin by creating 
a set of random positions with 5x the number of clusters 
as in the parent catalogue. We then calculate a clustering 
amplitude for each individual cluster in this random sample 
and attribute a weighting based on the assigned clustering 
amplitude. Clusters are then selected or discarded based on 
this weighting, until the realisation has the required number 
of selected clusters (i.e. that of the parent catalogue). As a fi- 
nal check, the auto-correlation of the realisation is measured. 
Through comparison with the auto-correlation of the origi- 
nal catalogue, the realisation is either accepted or discarded. 
This process is repeated until we have 100 acceptable clus- 
tered mock catalogues for each parent catalogue. The cross- 
correlation is then calculated between the WMAP data and 
the 100 mock catalogues and the standard deviation is taken 
as the lcr error on our results. 
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Figure 3. Cross-correlation results between 606 ACO rich galaxy clusters (R 2, |6| > 40deg) and the WMAP 3-year maps in 5 
band-passes (+ILC) as indicated. 



In addition to this we also perform a rotational analy- 
sis to provide an alternative estimate of the errors. In this 
case we perform the cross-correlation between the cluster 
positions and the WMAP data. We then shift the clus- 
ter positions by 20deg in galactic longitude and recalculate 
the cross-correlation. We repeat this until we have rotated 
through a full 360deg. A S/N is then calculated from the 
results of this rotational analysis. 



5 RESULTS 

5.1 Optical/IR Cluster Samples 

The results for the cross-correlation between the four large 
cluster datasets (APM m>7, APM m>15, ACO R>2 and 
2MASS clusters) and the WMAP W-band data are shown 
in Fig. [T] A decrement is immediately evident on small 
scales within 9 < 30' of cluster centres in all four data 
sets. Looking in detail first at the ACO results, the WMAP3 
cross- correlation strongly confirms the results of lMvers et al.l 
from WMAP 1st year data. Here, we find a decrement 
of -0.021 ± 0.007mK at 6 < 6'.3 and -0.010 ± 0.0Q4mK 
at 9 < 60' for the W-band data (quoted accuracies are from 
the Monte-Carlo analysis). Basically, the ACO decrement 
has remained the same and the improved statistics at small 
angles has increased the S/N. In addition to the Monte- 
Carlo analysis, we also checke d our ACO results using the 
rotational analysis described by |Mvers~e t al. (2004) and find 
the significance of the decrement at 6/3, 60' and 500' to be 
3.2a, 2.0cr and 1.2cr (see Fig. [2}. As before, there appears to 
be so me form of extended signal out to angles of ~ 100' . Fol- 
lowing [Mjers^etjil] l|2004h . we also produce the correlations 
with the four other WMAP bands, plus the ILC map. The 



results of this are shown in Fig. [3] Again, good agreement 
is seen between these updated results and the original 1st 
year data results. Despite the increasingly poor resolution 
of the bands, the decrement is observed in the V, Q and Ka 
bands, whilst even the ILC map and the Ka band map show 
a decrement. 

Improvements in the small scale statistics are also ob- 
served in the 2MASS and APM results while the magnitudes 
of the decrements remain unchanged. However, the APM 
group (m ^ 7) SZ detections remain marginal even at small 
scales. 



5.2 ROSAT X-ray bright cluster sample 

We next cons i der t he X-ray bright clusters of 
iBonamente et al] {2002). Analysis of this dataset with 
respect to the first year WMAP results has already been 
performed by iLieu et all (|2006l ). Their main conclusion 
was that the SZ decrement in the WMAP1 data around 
the locations of these clusters has a lower magnitude than 
they would expect from thei r predictions bas e d on the 
original X-ray observations of IBonamente et alj (|2002T ). In 
Fig. [4] (left panel) we show our c r oss-co rrelation between 
the 31 clusters used by ILieu et alj l|200rj) and the WMAP 
year 3 data in the W-band (crosses). W e also present the 
averag e model prediction based on the IBonamente et ahl 
(|2002h data (solid line). This has been convolved with a 
Gaussian beam pr ofile of a = 6/3. W e see the same general 
effect as seen by ILieu et all (|2006l ). that the SZ effect is 
somewhat smaller than predicted by the data. However, 
the significance of rejection is only w 2er. Similar results are 

seen in the other WM AP bands. 

We next split the lLieu et alj (|2006T > clusters by redshift 
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Figure 4. Average AT (from WMAP W-band data) plots for 
30 clusters from the ROSAT sample (left) and 39 clusters from 
the Chandra sample (right). In both figures, the points show our 
cross-correlation results, whilst the curves s how average SZ m od- 
els (based on the para meters taken from iLieu et al,l |2006| and 
iBonamente et al"1l2006h convolved with a Gaussian representing 
the WMAP beam profile. For the Chandra sample, we plot the 
full isothermal model (solid line) and the same model limited to 
9 < 2' (dashed line). 





Figure 5. Average AT (from WMAP W-band data) correlations 
for the ROSAT X-ray clusters split by redshift: 21 clusters at 
z < 0.1 (left) and 9 at z > 0.1 (right). 



as shown in Fig. [5] Again our results are given by the crosses, 
whilst the solid lines show the average model SZ profiles. 
The model for clusters at z < 0.1 is rejected by 4.2cr at 
8 < 6. '3, whilst at 0.1 < z < 0.3 the rejection drops to 1.6a 
at 6 < 6. '3. We have also performed this analysis with a 
latitude split at |6| = 40° and find some degeneracy between 
latitude and redshift as many of the low redshift clusters are 
also at low latitude. However, in either case we do not regard 
the difference between the results in Fig. [5] as particularly 
statistically significant. 



5.3 Chandra X-ray bright cluster sample 

We next analysed the S Z decrements for the 38 clusters of 
IBonamente et alJ |2006) , using the WMAP3 W band results. 
In Fig. U (right panel) we compare the cross-correlation re- 
sults with an average model constructed from the individ- 
ual is othermal models given in Table 5 of IBonamente et al.l 
l200d (solid line) and again find that the SZ effect is now 
quite severely over-predicted by the models, with a rejec- 
tion significance of 5.5a. We again looked for a dependence 
on redshift and found slight evidence for a greater SZ signal 
at z < 0.3 com pared to to z > 0.3 (Fig. El). 

Given that IBonamente et al.l l|2006h only fit the Chan- 
dra data for 6 < 2', there is the possibility that this model 
may not apply at the large angles covered by the WMAP 
data. We therefore also show in Fig. [4] the SZ model trun- 
cated at 8 = 2' before being convolved with the W-band 
beam (dashed line) . The significance of rejection in this case 





Figure 6. Top: Cross-correlations of WMAP Ka-band AT data 
with 20 clusters at z < 0.3 (left) and 19 at z > 0.3 (right) from 
the Chandra cluster sample. Bottom: The sa me for WMAP W- 
band data. The solid lines show the f3 models of Bonamente et all 
i2006l) convolved with the WMAP profiles. 



is reduced to 2.5cr. We note that this is a strict lower limit 
to this significance limit as it assumes no SZ contribution 
beyond 2'. 

Although within the Chandra sample there is little ev- 
idence of redshift dependence, the fits of SZ models to the 
WMAP data do appear to deteriorate as we move from the 
average redshift, z w 0.1 of the ROSAT sample to z ~ 0.3 of 
the Chandra sample (Fig. [4j| . We have also noted that at the 
lowest redshift the WMAP SZ effect is clearly detected at 
about the predicted amplitude in the Coma cluster (Fig. [7}. 
We therefore returned to the ACO dataset and identified 
407 R ^ 2, |6| > 40deg clusters with measured redshifts. 
Splitting these at z=0.15 (Fig. |8}, we see that there is some 
evidence confirming that clusters at higher redshift have ob- 
served SZ decrements that are significantly smaller than at 
lower redshift. Although the X-ray properties for the major- 
ity of these clusters are unknown, we have fitted the same 
average model, scaling 8 C to the appropriate average redshift 
before convolving with the WMAP beam. The fit appears 
significantly worse for the higher redshift clusters, with a 
rejection confidence at 8 < 6. '3 of ~ la for z < 0.15 and 
~ 4cr for z > 0.15. We tentatively conclude that there may 
be a redshift dependence of the SZ effect in the sense that 
higher redshift clusters show a smaller than expected effect. 



6 DISCUSSION 

The reduced SZ decrements in the WMAP3 data towards 
the ROSAT cluster sample and t he almost lack of detec- 
tion of the SZ effect in terms of the lBonamente et alJ (|2006l ) 
clusters is paradoxical. The most obvious explanation is that 
the WMAP data is contaminated by unresolved cluster radio 
sources within the WMAP beam. However, the contamina- 
tion from synchrotron radio point sources varies with fre- 
quency as T„ oc v a ~ 2 (where a ~ 0.7), whilst the d i screp - 
ancy in the WMAP3 data for the IBonamente et all (2006) 
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Figure 7. Binned AT data from the WMAP year-3 W-band 
data around the Coma cluster. The solid line shows the model 
predicted from X-ray data (taken from lLieu et al .120061) convolved 
with the W-band beam profile. 




Figure 8. Average AT (from WMAP W-band data) plots for 
the data from the Abell cluster catalogue with 172 clusters at 
z < 0.15 (left) and 235 at z > 0.15 (right). Only clusters with 
|b|^ 40deg are included here. Overlaid in both cases is a model 
with AT(0) = -0.16m.ftr, (3 = 0.7 and 8 C = 9/8 scaled to the 
mean redshift of the samples: z = 0.1 and z = 0.2. In both cases 
the model is convolved with the W-band beam profile. This model 
gives a reasonable fit to the data at z < 0.15, but significantly 
overestimates the z > 0.15 data. 



cluster sample is as large at Ka (33GHz) as at W (94GHz) 
(see Fig.©. 

Further to this issue, we note that a survey of radio 
sources in t he Ch andra clusters has been performed by 
ICoble et all (|2007l ). They see a population of radio sources 
with a mean flux of ~ 6. 2m Jy /cluster at 30GHz. Given a 
spectral index of such sources of a « 0.7 (where Sv oc v~ a ) 
this gives a flux of « 3m Jy /clust e r at t he W-band frequency 
of 90GHz. Following [Lieu et all (|2006t ). the equivalent flux 
required to cause the lack of SZ effect observed in the Chan- 
dra clusters can be determined from the Rayleigh- Jeans flux 
multiplied by the solid angle: 



Ssz = 



2nkATv 6 2 



(•>) 



Taking AT = O.lmK, v = 90GHz and 9 = 10', we ob- 
tain a flux of Ssz = 170m Jy. Even taking a value for 
the spectral index of a — for the radio sources (e.g. 
iBennett et al.l [ 2003i ), the flux required is over an order of 
magnitude greater than t he ob served discrete radio source 
fluxes from ICoble et all (|2007T l. In addition, iLin fc Mohr 
l|2007f ) make estimates of the contamination from radio point 
sources and for cluster masses typical of the Chandra sample 



(M200 ~ 1Q 15 M . iReiprich fc Bohringerl |2002|) . they sug- 
gest that up to only 10% of these clusters may be lost due 
to point source radio contamination (see their Fig. 15). Al- 
though this assumes that there will be no increase in source 
contaminatio n with the WMA P beam area, we note that 
the count s of ICoble et al.l (120071) are usually lower than pre- 
dicted bv lLin fc Mohrl 1 20071 . Fig. 13) and these effects may 
cancel. 

Currently we have no explanation for the strong SZ 
decrements detected by the interferometric experiments as 
opposed to the lack of detections by WMAP. At higher res- 
olution it may be more possible to detect the SZ against 
the noise caused by the primordial CMB fluctuations but 
our error analysis should take care of such statistical effects 
and the average model is rejected at the 5.5a level by the 
WMAP data. A high value of H « lOOkms'^pc -1 for the 
SZ X-ray model might help explain the ROSAT cluster re- 
sults but an even higher value would be required to explain 
the Chandra cluster results. 

As noted above there may also be evidence that the SZ 
de crement is too low in the ACO-WMAP1 cross-correlation 
of iMvers et~ai1 (|2004l ), as c onfirmed by t he AC O-WMAP3 
cross-correlation in Fig. [T] IMvers et al.l (|2004l ) noted that 
the decrement that fitted the ACO R ^ 2 clusters with = 
0.75 was only AT(0) = 0.083mK compared to the 0.5mK 
predicted for the R = 2 Coma cluster. The WMAP3 data 
confirms that AT(0) = 0.5mK is needed to fit the observed 
Coma SZ decrement (see Fig. [7}. In Fig. [T]the SZ models 
for these two values of the decrement are compared to the 
WMAP3 W band data for the R ^ 2 cluster sample. Both 
models assume j3 = 0.75. We see that while the data is 
well fitted at < 10' by the AT(0) = 0.083mK model, 
the AT(0) = 0.5mK at least begins to improve the fit at 
larger scales. One possibility is that, as well as detecting an 
extended SZ component to the ACO data, we may actually 
be detecting a lower central SZ amplitude than expected 
fro m the X-ray data. 

iLieu et ail l|200fj| ) discussed other possible explanations 
for the unexpectedly small SZ decrements det e cted in the 
ROSAT sample. For example, ILieu fc Quenbvl |2006) have 
discussed whether a diffuse cluster synchrotron source could 
explain the reduced SZ decrement. The main problem here 
is that non-thermal electrons would not give a good fit to 
the X-ray data which ar e usually well fitted by thermal 
bremsstrahlung, although ILieu fc Quenbvl (|2006l ) also noted 
that the soft X-ray excess seen in the central regions of some 
clusters may be indicative of a significant embedded non- 
ther mal X-ray componen t there. 

iFosalba et ak I (|2003l ) have discussed whether ISW effect 
could mask the SZ effect but the ISW effect is at 0.5/iK and 
seems too small to mask the SZ effect which in the X-ray 
clusters can be 10 x higher. 

There is also the possibility that the SZ decrement 
has been overestimated by the X-ray modelling. Certainly 
the Chandra predicted decrements for the 5 clusters in 
common with the ROSAT sample (A665, A1413, A1689, 
A1914, A2218) are on average w 80% larger than the pre- 
dicted decrements from the ROSAT data. Most of this dif- 
ference arises from A2218 where the ROSAT data imply 
AT(0) = -0.27mK (corrected to 30GHz) and the Chandra 
data imply AT(0) = -0.87mK, a factor 3.2 different. But 
since Chandra has higher spatial resolution, it is expected 
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to probe the central core of a cluster more accurately and 
so the Chandra X-ray models might be expected to be more 

robust than those from the ROSAT data. 

W hile this paper was in preparation, lAfshordi et al.l 
(2006) have also used X-ray data of 193 Abell Clusters to 
search for the SZ d ecrement from WMAP3 data (see also 
lAfshordi et al]|2005h . These authors made a significant de- 
tection and also suggested that the size of SZ decrement 
implied that the cluster hot gas fraction was 32 ± 10% 
lower than the baryon fraction in the standard cosmolog- 
ical model. They also suggested that their WMAP results 
were consistent with the interferometric SZ results for the 
sample of 38 C handra clusters analy sed above. Note that 
the approach of lAfshordi et all (|2006T l is different from that 
used here in that the X-ray data is mainly used to define a 
template to detect SZ decrements and then the SZ data and 
the X-ray temperature data alone are used to establish the 
gas densities. This route therefore avoids comparing the SZ 
results with X-ray gas density models on the grounds that 
the latter depend on assumptions such as that of hydro- 
static equilibrium. These authors also do not consider the 
possibility that the cluster SZ decrements may evolve with 
redshift. 

Finally, if we assume that the WMAP SZ decrements 
are reliable, even in the case of the 38 Chandra clus- 
ters where the unexpla ined discrepancy persist s with the 
OVRO/BIMA results of iBonamente et all l|2006h . we might 
speculate whether a lower than expected SZ decrement in 
the higher redshift clusters could be caused by foreground 
lensing. The indication from WMAP that the higher redshift 
clusters may have reduced SZ decrements is consistent with 
the idea that gravitational lensing is having a significant ef- 
fect on the detection of the SZ effect. Therefore it may be 
that the groups and clusters out to redshifts in the range 
0.2 < z < 0.8 in the foreground of the targeted Chandra 
clusters are lensing the cluster centres and smoothing the 
decrement away. Using CMBFAST we have constructed the 
lensing smoothing function for CMB scattering at z = 0.3 
and find that on the size of the » 10 arcmin WMAP beam, 
the smoothing function is reduced by about a factor of ~ 10 
compared to the case where the surface of last scattering 
is at z = 1100. At z — 0.7, the factor is ~ 5. Therefore 
for the standard model this would make the effect negligi- 
ble because at z=0.3, e = a/8 ~ 0.004 and at z — 0.7, 
e = a /6 m 0.008. Only if the mass power spectrum is sig- 
nificantly higher than that for the standard model can this 
explanation apply. One such case is the high mass power 
spectrum advocated by Shanks (2006) as a route to modify 
the first acoustic peak in the CMB. Such a spectrum is mo- 
tivated by the evidence from QSO lensing that the galaxy 
distribution is strongly anti-biased (b » 0.1) at least on 
0.1 — lh -1 Mpc scales with respect to the mass JMv ers et al.l 
l200ll2005l ; lMountrichas fc Shanksil2007h . However, the bal- 
ance of other evidence may still argue against such a high 
amplitude for the mass power-spectrum. 

Lensing would clearly also affect the X-ray cluster pro- 
files as well as the SZ decrements. Although these are ex- 
pected to be smoother than the SZ decrements, it might be 
expected that the profiles of lower redshift clusters are on 
average steeper than the profiles of higher redshift clusters. 
It remains to be seen whether this prediction of the lensing 
hypothesis can be decoupled from evolution of the cluster 



gas component. In any case, the flatness of the X-ray pro- 
files towards the centres of many clusters may make this 
prediction more difficult to test. 



7 CONCLUSIONS 

We have confirmed the extended appearance of the SZ decre- 
ment in WMAP 3-year da ta around ACQ R ^ 2 clusters out 
to 9 » 30', first shown bv lMvers et~ai1 l|2004h using WMAP 
1-year data. Further to this, we have confirmed the detec- 
tion of the SZ decrement in the 3-year data around clusters 
identified in both the APM survey and 2MASS, showing 
an increase in detection significance compared to the 1-year 
data analysis. 

We have also confirmed the result of iLieu et al.l (2006) 
that the SZ decrement is somewhat lower than expected 
on standard model assumptions and R OSAT X-ray pro- 
files f or a sample of 31 clusters from IBonamente et all 
(2002). We have further shown that even smaller X-ray 
decrements are seen in the higher -redshift sample of 38 
clusters of IBonamente et all (|2006l ) that has Chandra X- 
ray data. The reason for the observational discrepancy be- 
tween the WMAP dat a and the BIMA/OVRO data of 
IBonamente et all (|2006T > is not clear. We do not believe that 
discrete or diffuse cluster radio sources nor the ISW effect is 
likely to explain the discrepancy. Dividing the ACO clusters 
into high and low redshift samples also indicates that the 
deficit in the SZ decrement may increase at higher redshift. 

In the light of the above results from our WMAP SZ 
analysis, we have discussed the possibility that the extended 
SZ signal detected for ACO and 2MASS clusters may actu- 
ally be indicating a lack of SZ signal in the centres of clusters 
rather than an excess at the edges. 

On the assumption that the WMAP SZ results are cor- 
rect, one explanation we have considered is that lensing of 
the cluster centres by foreground groups and clusters could 
explain the over-prediction of the observed decrements by SZ 
models and in particular the apparent tendency for higher 
redshift clusters to have smaller SZ decrements. However, 
before considering such interpretations further, we need to 
clarify if this is a real observational discrepancy between the 
OVRO/BIMA data and WMAP. 

It will clearly be interesting to see if these WMAP re- 
sults are confirmed in the higher resolution SZ observations 
made using the Planck satellite. 
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